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Abstract 



We study an extended Lipkin-Meshkov-Glick model that permits a transi- 
tion to a deformed phase with a broken continuous symmetry. Unlike simpler 
models, one sees a persistent zero-frequency Goldstone mode past the tran- 
sition point into the deformed phase. We found that the RPA formula for 
the correlation energy provides a useful correction to the Hartree-Fock en- 
ergy when the number of particle N satisfies > 3, and becomes accurate 
for large N. We conclude that the RPA correlation energy formula offers a 
promising way to improve the Hartree-Fock energy in a systematic theory of 
nuclear binding energies. 
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I. INTRODUCTION 



Hartree-Fock (HF) theory is the fundamental starting point to understand the ground 
state properties of many-fermion systems. Its main assumption is that a particle indepen- 
dently moves in a mean field generated by other surrounding particles. In nuclear physics, 
with the adjustment of a few parameters of an effective interaction, the HF theory has 
described reasonably well the global properties of nuclei throughout the periodic table [|1|. 
However, correlation effects which go beyond HF are also significant. They appear most 
dramatically when the HF ground state violates a symmetry of the Hamiltonian such as ro- 
tational invariance or number conservation. A global theory of nuclear binding must surely 
take these correlation into account, if it is to achieve an accuracy at the 1 MeV level. A 
recipe is often used which is based on the projection after variation method [0. 

In this paper, we argue that the RPA correlation formula can provide a better method 
in correcting the broken symmetry in the HF theory. In the RPA, a restoration of the 
symmetry breaking appears as a zero energy solution of the RPA equations. The correlation 
energy associated with the several symmetries can be calculated using the RPA formula 

Ecorr = ^ (E - Tr{A)^ , (1) 

where Ui is the (positive) frequency of the RPA phonon for the z-th mode and A is the A 
matrix in the RPA equations. Our aim in this paper is to construct a simple Hamiltonian 
model to study the effects of correlations on the ground state energy and to show that the 
RPA formula is adequate in principle. To this end, we employ a three-level version of 
the Lipkin-Meshkov-Glick (LMG) model. The model describes identical fermions in three 
single-particle levels, each of which is A^-fold degenerate. Exact solutions can be obtained by 
explicitly diagonalising the model Hamiltonian. The RPA correlation energy was discussed 
by Brito et al. 0, but their parameters did not leave a continuous symmetry to be broken, as 
is the case for the transition between spherical and deformed nuclei. In this paper, we shall 
use parameters which mimic quadrupole motions in nuclei. In such a way, the correlation 
associated with rotational motion is easily studied. 
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The paper is organised as follows. In Sec. II, we first show that the RPA formula is 
exact for a simple model with a two-body interaction, Eq. (|^) below. We then introduce 
the three-level Lipkin model in Sec. Ill and solve it in the HF as well as in the RPA. We 
compare the RPA correlation energy with the exact solution of the model and show that 
the RPA formula works well even in the vicinity of the critical point of the phase transition. 
Summary of the paper is given in Sec. IV, together with further discussions on the RPA 
formula. 



II. RPA CORRELATION FORMULA 

Before we study the correlation in the three-level Lipkin model, we would like to demon- 
strate that the RPA correlation formula works well, using an analytically solvable model. 
Consider a two fermion system bound in a harmonic potential coupled by a linear interaction. 

This model was first introduced in Ref. to discuss the effects of the RPA correlation on 
the ground state density. A similar model has been considered in Ref. P] in connection with 
the paired Wigner crystal. Using the transformation ^ = (xi +X2) / \/2 and rj = (xi — X2) / V^, 
the Hamiltonian can be written in the form 

from which we obtain the exact ground state energy 



Eexact = lhuJo{Jl- + Jl + -^] . (4) 

The limit where the coupling constant C is equal to mt^Q corresponds to a translationally 
invariant Hamiltonian with an interaction ^^mujQ^Xi — ^2)^, giving 

Eexact = Ti^o/ V2- (5) 
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Let us now solve the problem in the mean field approximation and then consider the 
RPA correlation energy. If the coupling constant C is small, one can regard V in the 
Hamiltonian as a residual interaction. The mean field Hamiltonian Hq then has already 
been decoupled and we immediately obtain 



Emf = fiuJo/2 + huJo/2 = huQ. 
We define the RPA excitation operator as 



(6) 



j=l,2 



(7) 



where and a are the creation and the annihilation operators of the unperturbative phonon, 
respectively. The RPA equation then reads 
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(8) 



ao being the amplitude of the zero point motion defined as ^Jh/2mujQ. The solutions of this 
equation are found to be 



TlUJ = ±1 TlOJn 
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We thus obtain 
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which is precisely the needed correction to get the ground state energy Eq. starting from 
the mean field energy Eq. (^). It reproduces Eq. in the translationally invariant case. 

It is interesting to compare the RPA approach with other ways of dealing with correlation 
energies associated with broken symmetries. In the case of center-of-mass motion, a recipe 
is often to subtract the expectation value of the center of mass operator from the mean field 
energy. With our Hamiltonian, this prescription gives 
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E,m = -<MF\-^^^^^\MF>=--huJo. (11) 

The total Emf + = SHluq/A is not exact, although it is close to Eq. (|^). This study 
clearly shows that the RPA formula provides a much better method to calculate correlation 
energies. 



III. THREE-LEVEL LIPKIN MODEL 

We now consider RPA correlations in a three-level Lipkin model. Labelling the levels 
0,1, and 2, we choose the Hamiltonian to be invariant under transformations between 1 and 
2. The Hamiltonian we consider can be expressed 

H = e(ni + ns) - jiK^K, + K2K2 + kIkI + kIkI), (12) 

where 

N 



na = ^C^aiCai, « = 0, 1, 2, (13) 

i=l 

N 

Kc, = Y.ciiCoi, a = 1,2. (14) 



i=l 



A. Exact Solutions 



Since the Hamiltonian given by Eq. (12) couples symmetric states with respect to inter- 
change of particles only with other symmetric states, a suitable basis for the exact diago- 
nalisation of the Hamiltonian H is given by [H] 



This is a simultaneous eigenstate of the number operators fii and 712 with the eigenvalue of 
Hi and 712, respectively. The effect of the K operators on the states is given by relations 
such as 



5 



Ki\nin2 >= \J{N - rii- n2){ni + 1)|?t,i + 1, ?t,2 > 



(16) 



The matrix elements of H can easily be calculated and are given by 



< n\n'2\H\nin2 > 



e{ni +Ti2)5n;,ni5n: 

V 



,n2 



-— (^\/ (ni + l){ni + 2){N - ni ~ n2)iN - rii - n2 - l)5„'^,ni+25ni,,n2 
+ \Jniini - 1){N - ni - 722 + 1)(A^ - ni - n2 + 2)5„'^,„^_2(5^,n2 
+ \l {n2 + l)(n2 + 2){N ~ni- n2){N - ni - n2 - l)5„4,ni5n^,n2+2 



+ \/«2(^2 - l)(iV - - n2 + 1)(A^ - rii - n2 + 2)5„'^,„^5„^,„2_2 

The dimension of the matrix to be diagonalised is (A^ + 1)(A^ + 2)/2. Further reduction can 
be achieved by considering that the Hamiltonian conserves the parity of each levels [§. 



B. Hartree-Fock Approximation 

Let us now solve the problem in the Hartree-Fock approximation. We consider a trans- 
formation of basis defined by operators a^i, with agj representing the occupied orbital. The 
HF state has the form 



N 



\HF>=Ylal\>. 



(18) 



and the transformation of basis is such as to minimize the expectation of the Hamiltonian. 
Let us write the transformation as 



( 



\ 0-21 ) 



\ 



COS a cos /5 sin a sin /5 sin a 
— sin a cos /? cos a sin fi cos a 
— sin /3 cos /3 

Using these relations, it is straightforward to evaluate the energy surface E{a,P) 
= < HF\H\HF > as 



V 



J 



\C2,J 



(19) 



E{a, p) = Ne sin^ a - VN{N - 1) sin^ a cos^ a. 



(20) 
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Note that the potential surface E{a, (3) is independent of (3 and thus totally flat in the 13 
direction for the rotationally invariant Hamiltonian. For simplicity, we particularly choose 
/? = in constructing the HF single particle operators, Eq. (|19|). The HF Hamiltonian thus 
spontaneously breaks the rotational symmetry, and the Goldstone mode will appear at zero 
excitation energy to restore the symmetry breaking, as we will show in the next subsection. 

The optimum choice of a is obtained by minimising the potential surface E{a,(3). It is 
convenient to express the solution in term of the dimensionless parameter 

X^V{N-l)/e. (21) 

For X < 1) the minimum appears at a = (spherical phase). At x = 1, the system undergoes 
a phase transition and, for x > the potential surface displays two symmetrical minima at 
cos2q; = 1/x (deformed phase). The ground state energy in the HF approximation is thus 
given by 

'o (x<l), 
Ehf = { (22) 

^(2-x-l/x) (X>1). 



C. Random Phase Approximation 

We next solve the problem in the RPA in order to evaluate the correlation energy asso- 
ciated with the rotational motion. We define the RPA excitation operator as 

= XiKi + X2K2 - YiKl - Y2KI (23) 

where 

_ N 

Ka = Y^ aLaoi, a = 1, 2. (24) 

i=l 

The RPA equation is obtained from < HF\[5Q,[H,Q^] - uQ^]\HF >= for 6Q = 
Ki , K2 , K\ , and K2 ■ The result is the well-known RPA matrix equation 
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(25) 



where A and B are 2x2 matrices given by 



On 

Au — e cos 2a + -ex sin 2a, 



Ai2 = A 



21 



0, 



A22 — e(l — sin^ a) + -e^sin^ 2a, 
Bii = — ex(cos'^ a + sin'' a), 
B\2 — B21 — 0, 



B. 



22 



-ex cos a. 



(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



Because A and B are separately diagonal, the RPA matrix can be easily diagonalized. The 
RPA frequencies are found to be 



{An + B^,){A^,-B^,), 



^2 — (^22 + -B22)(^22 — -622)- 



(32) 
(33) 



Substituting the self-consistent value for a obtained in the previous sub-section, we obtain 



(34) 



(35) 



e2(l + x)(l-x) (X<1), 

2e'(x + l)(x-l) (X> 1), 

e2(l + x)(l-x) (X<1), 

(X > 1)- 

In the spherical phase, the RPA frequencies for the two modes are identical. In the deformed 
phase, on the other hand, the frequency for the second mode becomes zero. In this case, 
the first mode corresponds to the beta vibration, while the second mode corresponds to the 
rotational motion perpendicular to the symmetry axis. Figure 1 shows the RPA frequencies 
as a function of x for N = 20. One can clearly see the discontinuity at the critical point 



Figure 2 compares the ground state energy as a function of x obtained by several methods. 
The number of particle is set to be 20. The solid line is the exact solution obtained by 
numerically diagonalising the Hamiltonian. The dashed line is the ground state energy in 
the Hartree-Fock approximation given by Eq. (p2D. It considerably deviates from the exact 
solution through the entire range of x shown in the figure. The dot-dashed line takes into 
account the RPA correlation energy in addition to the HF energy. The improvement is 
apparent and significant. It is remarkable that the RPA formula works well even in the 
vicinity of the critical point of the phase transition x = 1- 

It is expected that the RPA should be accurate for a large number of particles, and 
A^ = 20 seems to fulfil that condition. We also would like to apply the approximation to 
nuclei with just a few valence nucleons. Figure 3 compares the exact energy with the HF 
+ RPA as a function of A^, for x=5.0. The RPA correction may be considered useful if it 
is within a factor of 2 of the exact. We see from the figure that this is satisfied for A^ > 4. 
For N = 2, the exact correlation energy is several times RPA, and such near-magic nuclei 
would require a more elaborate way. 

IV. SUMMARY AND DISCUSSIONS 

We discussed the role played by the RPA correlation in the ground state energy. To this 
end, we used simple Hamiltonian models, like a bi-linear interaction between two fermions 
as well as a three-level Lipkin model. The former is for the correlation associated with the 
center of mass motion, while the latter for the rotational motion. We showed that the ground 
state energy is well described in the mean field theory once the ground state correlations 
are taken into account in the RPA. We also showed that the RPA formula works well for a 
wide region of an order parameter, including in the vicinity of the critical point of the phase 
transition. Evidently, the RPA formula provides a powerful method to calculate energies for 
long range correlations, which are not included in the HF approximation. 

Up to now, microscopic theory based on the mean field theory has not been as successful 
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than other approaches in making a global fit to nuclear binding energies. The most accurate 
theory of nuclear binding systematics |Ty] starts from the liquid drop model, and treats shell 



effects perturbatively. It fits the binding energies with an RMS deviation of 0.67 MeV, a 
factor of three better than the Gogny or the published Skyrme functional. We believe that 
there are good prospects to develop a better microscopic global theory, treating correlation 
energies systematically by the RPA. A work towards this direction is now in progress [[Till . 
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Figure Captions 

Fig.l: RPA frequencies for collective vibrations as a function of x = y(A^ — l)/e. The 
number of particle N is chosen to be 20. 

Fig.2: Comparison of the ground state energy obtained by several methods. The solid hne is 
the exact numerical solution. The ground state energy in the Hartree-Fock approximation is 
denoted by the dashed line, while the dot-dashed line takes into account the RPA correlation 
energy in addition to that. 

Fig.3: The ground state energy as a function of N ior x — 5.0. The meaning of each hne 
is the same as in Fig. 2. 
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